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Avalanche photodiodes are commonly used as detectors for low energy x-rays. In this work we report on 
a fitting technique used to account for different detector responses resulting from photo absorption in the 
various APD layers. The use of this technique results in an improvement of the energy resolution at 8.2 keV 
by up to a factor of 2, and corrects the timing information by up to 25 ns to account for space dependent 
electron drift time. In addition, this waveform analysis is used for particle identification, e.g. to distinguish 
between x-rays and MeV electrons in our experiment. 


I. INTRODUCTION 

Avalanche photodiodes (APDs) are silicon-based solid 
state detectors that convert photons into a charge cur¬ 
rent. They provide a compact, robust, magnetic field 
insensitive solution for light and x-ray detection with 
gains on the order of 100 and fast response times®^^. Due 
to this, APDs are extensively used in a large variety of 
physics®^®, medical^ and aerospace applications®®. 

We have studied x-rays with energies between 1 and 
lOkeV and observed two distinct APD responses to mo- 
noenergetic x-rays absorbed in different depths inside the 
APD. By constructing APD specific standard traces, and 
using a pulse-by-pulse fitting technique, we improved the 
APD energy resolution by a factor of 2, and the time 
resolution by 30%. In addition, we were able to identify 
background signals stemming from electrons that deposit 
a few keV energy in the APD. 

The data presented in this work were gathered in the 
muonic helium Lamb shift experiment®®’®^, using a set 
of twenty large area avalanche photo diodes (LAAPDs) 
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from Radiation Monitoring Devices (model S1315; 13.5 x 
13.5 mm® active surface area each). The muonic he¬ 
lium ions represent an extended x-ray source that emits 
predominantly monoenergetic x-rays of 1.52 keV and 
8.22 keV as well as electrons with up to 50 MeV of kinetic 
energy (see Appendix). Previous tests of these APDs 
found 40% detection efficiency for 8.2 keV x-rays, and an 
average energy resolution of 16% (FWHM) after calibra¬ 
tion®®. 

Our x-ray detection setup consists of two linear arrays 
of 10 LAAPDs each, in which each LAAPD is mounted 
on a separate titanium piece for efficient cooling and easy 
replacement®’®®. The detector arrays are mounted inside 
a vacuum around 10“® hPa, and inside a 5 Tesla mag¬ 
netic field, above and below the x-ray source. Custom- 
built low-noise, fast response preamplifiers are fitted to 
the LAAPDs. Both LAAPD/pre-amplifier assemblies are 
cooled using an external ethanol circulation system and 
are actively temperature stabilized at around —30°C. 
The achieved short term temperature stability was bet¬ 
ter than ±0.1°C. Highly stable temperatures are crucial 
for the operation of LAAPDs since their gain depends 
strongly on their operating temperature®’®®. Bias volt¬ 
ages were chosen to provide the best energy resolution 
per APD and ranged from 1.61 kV up to 1.69 kV, approx¬ 
imately 50 Volts below the breakdown voltage. The pre¬ 
amplifiers with two bipolar input transistors in cascode 
configuration (BFR 182 npn, BFT 92 pnp) have been 
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FIG. 1. X-ray energy spectrum from a single APD before 
(top) and after (bottom) applying onr correction. The first 
spectrum is obtained by integrating over the recorded pulse 
amplitude in a 200 ns time window after the leading edge. A 
difference in extracted energy for detector responses with slow 
and fast rise times (labeled slow 8 keV and fast 8 keV x- 
rays respectively) is clearly visible. Improved energy calibra¬ 
tion managed to unite both responses and improve the energy 
resolution by up to a factor of 2 (from 32% to 16% FWHM at 
8.2keV for this APD). 


used for the generation of a fast response from the large 
capacitance (120 pF) of the LAAPD. An overall gain of 
150 mV /fjiA at 5011 has been measured with a test pulse. 
Outgoing APD signals were further amplified by gain 4 
main-amplifiers and fed to the CAEN vl720 waveform 
digitizers (250 MS/s, 12 bit) for recording. 

Our experiment requires pileup detection in the x-ray 
detectors to reduce background effects. Standard shap¬ 
ing amplifiers that are commonly used feature integration 
times too long to separate pulses on a 100 ns scale. This 
deteriorated the performance in our previous measure- 
ments^^“^® where we used Rutherford Appleton Labora¬ 
tory (RAL) 108A pre-amplifiers with /rs-long integration 
times (see^^ Fig. 16). For our new project^^d 2 ^gg^j 
fast pre-amplifiers with 30 ns rise time. When calculat¬ 
ing a simple integral over the recorded pulses, a poor 
energy resolution became visible as seen in Fig. 1. The 
double peak structure that was clearly resolved in 6 out 
of 20 APDs is a result of two different APD responses to 
the monoenergetic 8.2 keV x-rays as can be seen in the 
upper part of Fig. 2. Similar effects were previously re¬ 
ported for beveled edge APDs and 14.4keV x-rays^^. We 
first observed the same behaviour in a separate test setup 
without magnetic field. Hence the features described here 
can not be attributed to magnetic trapping effects in the 
drift region^®. We can only speculate why this effect was 
not seen for x-rays taken in another experiment at cryo¬ 
genic temperatures^^. Pre-selection of APDs with good 



FIG. 2. Top: Typical APD responses for 8.2 keV x-rays. Even 
though incoming x-rays are quasi mono energetic, the APDs 
show two distinct responses. The fast 8 keV component has 
a rise time of about 35 ns while the slow 8keV component 
shows a rise time of about 70 ns. Separate averaging of both 
individual data sets allows to produce standard traces that 
accurately describe all x-ray traces between 2 keV and 10 keV. 
Bottom: Electron induced signals that correspond to an x-ray 
energy of 8.2 keV after calibration. The dashed curves show 
the average of the slow 8 keV and fast 8 keV x-rays. Even 
though similar in shape to fast 8keV x-ray signals, a fit 
was able to identify 86% of these electrons correctly. 


energy resolution at 5.9 keV can lead to the vanishing 
of the double peak structure. This was the case in our 
previous measurement^®. Also the large average angle of 
incidence in our setup increseases this effect significantly. 

To compensate, we developed a simple standard re¬ 
sponse fitting technique that allowed us to distinguish 
between different responses on a hit-by-hit basis, improv¬ 
ing the energy resolution by a factor of two (see Fig. 1, 
bottom) and correcting for a 25 ns time shift between 
both signal types as discussed in Section IV. 

In the Secs. II-VI, the different features of the mea¬ 
sured x-ray signals are discussed before the fitting rou¬ 
tine and the improved energy calibration are presented. 
Then timing difference between both responses and the 
influence of electron signals in the analysis are reviewed 
before a brief summary and outlook is given. 


II. APD X-RAY RESPONSE 

The working principle of the APDs used in our setup is 
explained in Fig. 3. In the conversion region (II), incom¬ 
ing photons produce primary photoelectrons. Differences 
in the thickness of this layer (II) give rise to changes in 
detector energy acceptance. A p-n junction is placed on 
the back side of the active volume creating high local 
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FIG. 3. Working principle of avalanche photodiodes based 
on a p’'"-i-p-n^ doping profile. The weakly doped intrin¬ 
sic part (II) serves as conversion region for most incoming 
x-rays (case 1). Photoelectrons created are transferred to¬ 
wards the avalanche region. In this high field area secondary 
electrons are generated through impact ionization providing 
charge gain. Low energy x-rays have a high probability of 
being stopped in the initial drift region (I) (case 2). These 
experience additional signal delay and reduced gain. Some 
photons convert in the multiplication region (III), also lead¬ 
ing to reduced signal amplitudes (case 3). More about this 
effect can be found in^^. The bottom figure shows the electric 
field profile in the several regions of the APD together with 
the x-ray absorption profile for 1.5 keV and 8.2 keV x-rays. 


field strengths. Inside this avalanche region (III) elec¬ 
tron impact ionization at the high field p-n”*^ junction 
leads to a multiplication of free charge carriers providing 
gain for the initially converted primary photoelectrons 
The calculated absorption length for 8.2 keV and 1.5 keV 
x-rays is 70^m and 9/xm, respectively^®. Due to the 
extended size of our x-ray source, the average incident 
angle of 52 degrees in our geometry gives rise to an effec¬ 
tive 1.6 times longer average path inside the APDs. This 
absorption length for 8 keV x-rays is similar to the APD 
layer thicknesses and therefore leads to a number of dif¬ 
ferent effects on the APD output depending on the region 
where the photon is absorbed. The different possibilities 
are also shown and explained in Fig. 3. 

The largest part of the recorded 8 keV x-rays stops in 
the conversion region (II) and follows the normal APD 
working principle that provides high charge collection ef¬ 
ficiency and fast amplification. Nevertheless, some x- 
rays are absorbed either in the drift layer (I) or in the 



FIG. 4. Normalized slope of the rising edge plotted versus the 
integral of the pulse. The z-axis (color scale) is logarithmic. 
Integrals are roughly proportional to the deposited energy of 
the registered x-rays. Four contributions are visible: Low 
energy 1.5 keV x-rays show integrals below 200. The recorded 
8.2 keV x-rays create two different responses in the APD, one 
with slow rise time (slope « 0.3), and one with significantly 
faster rise time (slope « 0.7). The last contribution with an 
integral above 700 arises from MeV electrons depositing keV 
energy in the APD active region. 


avalanche region (III). The x-rays absorbed in region 
(III) undergo only partial amplification resulting in low 
amplitudes down to zero. This gain reduction is respon¬ 
sible for the fiat energy tails seen in Fig. 1. X-rays ab¬ 
sorbed in region (I) generate electrons which are only 
slowly transfered to the following region (II) due to the 
lower field strengths in (I). Traps in this region may 
hold electrons for non-negligible times, lengthening the 
pulse and causing a reduction in amplitude (see Fig. 2). 
Similar effects of reduced charge collection efficiency were 
also studied for x-ray energies below the silicon K-edge 
From Fig. 2 we also observe that these x-rays only show 
a single amplitude and not a continous distribution up to 
the one of x-rays absorbed in region (II). This indicates 
that the trapping mechanism occurs at the boundary be¬ 
tween regions (I) and (II). 

III. X-RAY ENERGY DIFFERENCES AND 
COMPENSATION 

In order to investigate this effect a set of roughly 
2.5 X 10"*^ x-ray traces were recorded per APD. Fitted 
baseline fluctuations were below 10 mV for all analyzed 
signals, compared to average signal amplitudes of 500 mV 
for 8.2 keV x-rays. 

Our analysis routine starts with an edge finder (square 
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FIG. 5. Energy spectra of recorded x-rays and electrons in 
the muonic helium Lamb shift experiment categorized by the 
standard trace which provides the lowest in n range of 
200 ns after the leading edge of the pulse. All spectra show 
two prominent peaks at 1.5 keV and 8.2 keV. The fast rising 
component provided by Trfast in dark blue enfolds signals 
converted in or behind the conversion region (II). It con¬ 
sists mostly out of 8.2 keV x-rays and the visible low energy 
tail is created by the loss of gain for x-rays converted in the 
avalanche region (III). The light blue distribution stands for 
all traces that were best described by the slow rising pulse 
shape TtsIow and consists mostly out of 1.5 keV x-rays and 
some 8.2 keV x-rays mixed in. Signals best matching the elec¬ 
tron trace TVeiec are shown in the orange division. These 
signals are formed by a continuous electron background and 
a contribution of wrongly identified x-rays. 


weighting function with a width of 200 ns) to find the 
beginning of the pulse in the recorded trace. Then the 
slope of the leading edge is fitted with a linear function. 
Using a criterion we improve the accuracy of the slope 
determination by varying start time of the pulse within 
20 ns while keeping the fitting window fixed. Finally, we 
normalize the slope to the pulse integral provided by the 
edge finder to obtain the (amplitude-independent) rise 
time of the pulse. 

When the rise time is plotted versus the integral of 
the pulse, four different contributions to the spectra can 
be identified as seen in Fig. 4. The two most promi¬ 
nent peaks are created by converted 8.2 keV photons with 
slow and fast detector responses, labeled slow 8 keV and 
fast 8keV respectively. For these peaks we see a clear 
difference in rise time and integral while most of the low 
energy 1.5 keV x-rays show a slow rise time. The rise 
time distribution for small signals is broadened due to 
low amplitudes and noise. 

The last visible component is generated by the al¬ 
ready mentioned high energy (up to 50MeV) electrons 
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FIG. 6. Spectra showing the relative timing between 
slow 8keV and fast 8keV x-rays for a single APD. These 
time spectra have been optained by plotting the time of the 
8.2 keV x-ray signal detected in a LAAPD relative to the 
1.5 keV signal detected in another LAAPD for both classes. 
Both LAAPDs including preamp-delay line etc. have been 
synchronized using electrons. The origin is chosen as the cen¬ 
ter of gravity of the slow 8 keV response. A 25 ns timing dif¬ 
ference could be measured between both responses. The rela¬ 
tively poor time resolution is given by the coincident 1.5 keV 
x-rays that are detected just above the noise level in our setup 
leading to a washed out signal. 


(created by muon-decay, further explained in the ap¬ 
pendix). These electrons deposit energies up to 50keV 
in the APDs and their signals display a third kind of 
standard pulse shape, namely a mixture of fast and slow 
x-ray pulse shapes. This is shown in the lower panel of 
Fig. 2. 

In order to further analyze the two classes of 8.2 keV 
x-rays, two sets of APD traces for slow 8keV and 
fast 8 keV were created by selecting the respective peaks 
in Fig. 4 with adequate cuts. A collection of selected 
traces for the slow 8 keV and fast 8 keV cases is shown 
in Fig. 2. For each of the two x-ray classes, traces were 
numerically averaged after shifting each trace to correct 
for the variation of the pulse starting time. This aver¬ 
aging created the standard traces of the subsets (TVsiow 
& Trfast)- These traces had to be produced once per 
each APD for a measurement period of several months 
and stayed constant throughout multiple heating/cooling 
cycles of the APD assembly. 

For the final analysis, each APD pulse is fitted with all 
available standard traces. Starting at the time provided 
by the edge finder, the standard trace is fitted to the 
pulse. The timing is then varied and the is recorded 
for each fit. To save computational effort that would 
arise for a 2-parameter fit (amplitude and time), the am¬ 
plitude of the standard trace is always fixed by matching 
its integral to the integral of the signal (after baseline 
subtraction) in a 200 ns wide time window. 

Finally, the minimal between the various standard 
traces is used to separate the pulses into different classes: 
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slow 8keV and fast 8keV (and electrons, see below). 
The result from the best-fitting class is used to get am¬ 
plitude, integral and timing values of the recorded signal. 

The allocation of the recorded APD signals in the 
slow 8 keV and fast 8 keV classes according to the fit 
routine can be seen in the top parts of Fig. 5. Calibration 
of the two x-ray spectra created by the TVsiow and TVfast 
fits is done by matching the peaks in both separate in¬ 
tegral spectra to the respective energy of 8.2 keV. As ex¬ 
pected, the fast 8 keV component of x-rays is the largest 
part of the recorded signals in our setup as seen in Fig. 5. 
The observed 1:1.7 ratio of fast 8keV to slow 8keV 
x-rays agrees roughly with the expected absorption ratio 
of 1:1.5 estimated from the thicknesses of layers (I) and 
(II). 

IV. X-RAY TIMING DIFFERENCES 

In addition to the variation in the observed 8.2 keV 
x-ray energy we were also able to measure a difference 
in timing between the fast 8 keV and the slow 8 keV 
components. 

In order to achieve a common timing reference point 
for this study, coincidence events between the 8keV x- 
rays recorded in the APD under investigation and the 
1.5 keV x-rays registered in neighboring APDs were stud¬ 
ied. These two x-ray types are emitted within a picosec¬ 
ond time window from the muonic atoms used (as is fur¬ 
ther explained in the appendix). 

Special attention was given to time calibration in or¬ 
der to avoid possible timing shifts created by the distinct 
standard traces TVfast and TVsiow used for different sig¬ 
nals. Therefore calibration of the different APDs and 
traces against each other was done using the supplemen¬ 
tary measured electron signals. The MeV electrons create 
hits in multiple detectors on their spiraling motion in the 
surrounding magnetic field enabling us to get a common 
timing for all APDs. 

When comparing the timing of the measured 8.2 keV 
x-rays we observed a 25 ns delay between slow 8keV 
signals and normal fast 8 keV signals. A time spectrum 
showing this effect for a single APD is shown in Fig. 6. 
Correcting for this effect improves the APD time res¬ 
olution of our setup by more than 30 % when the two 
responses seen in Fig. 6 are unified. Better results might 
be achieved when a more clearly defined common timing 
is provided since the timing resolution is limited by the 
low amplitude 1.5 keV signals just above the noise level. 

V. MeV ELECTRON DETECTION WITH APDS 

Apart from the improved energy resolution that was 
achieved with the methods described in Section III, we 
were also able to differentiate high energy electron sig¬ 
nals in the APDs from similar x-ray signals. These MeV 
electrons deposit up to 50 keV in the APD active volume 


and were always present in the experiment. Due to their 
passage through all the APD layers, electrons show sig¬ 
nals with yet another shape that can be distinguished 
from the previously discussed fast and slow x-ray re¬ 
sponses. A third standard trace TVeiec was created by 
averaging a set of clearly identified electron signals that 
correspond to a mean energy of 12 keV. This was done 
supplementary to the already known x-ray traces TVsiow 
and TVfast- A comparison of an electron induced signal 
shape at 8.2 keV and the respective real x-ray traces is 
also shown in Fig. 2 (bottom). As electrons with MeV 
energies deposit energy in all three APD layers. (I, II, 
III). The corresponding standard trace TVeiec can be 
approximately parameterized as a mixture of the TVgiow 
and TVfast standard responses. Using the same routine 
as for the previous pulse analysis, the fit was able to dif¬ 
ferentiate between x-ray and electron signals with very 
high fidelity, leading to a correct electron identification 
in 86% of the cases. 


VI. SUMMARY 

We have observed effects from the APD layer struc¬ 
ture that lead to two distinct responses to X-rays in the 
6-10 keV range. The individual signal types can be iden¬ 
tified with high fidelity by examining the rising edge of 
the measured pulses. Correcting for this effect improves 
the energy resolution by up to a factor of 2 depending 
on the APD. Additionally we were able to correct for 
timing differences between both responses. While the 
different rise time classes were observed in all 20 APDs 
under investigation, only 6 of them showed a resolved 
double-peak structure in the energy spectrum obtained 
by a simple integral. 

Using the rise time analysis, it was also possible to filter 
MeV energy decay electrons. An electron-specific stan¬ 
dard trace was clearly distinguishable from the two dif¬ 
ferent kinds of x-ray signals recorded for 8.2 keV x-rays. 
A fit of the signal shape was used to exclude them 
from the x-ray data with an overall effectiveness of 86 %, 
while only 14 % of the 8 keV x-rays were wrongly iden¬ 
tified as electrons. This lead to significant background 
reduction in the /xHe Lamb shift experiment^^’^^. 
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Appendix: The /rHe Lamb shift experiment 

The data presented in this work were acquired using 
muonic helium ions as x-ray source during the recent 
/iHe Lamb shift experiment^^d2 experiment is per¬ 
formed at the high intensity proton accelerator facility 
at Paul Scherrer Institute in Switzerland. Its’ purpose 
was to measure the different 2S —> 2P transitions in the 
/i^He"'' and /i^He'*’ exotic ions via laser spectroscopy. The 
required information about its environment and working 
principle will be briefly sketched in this section. 

The accelerator physics environment leads to stringent 
demands on stability and robustness of the APDs and 
the analysis routine employed that exceed common spec- 
ihcations. For example, the APD arrays used are placed 
inside a 5 T solenoidal magnet where they are mounted 
next to a low pressure helium gas target. Muonic ions 
are created in this 20 cm long gas volume operated at 
2 — 4 hPa by low energy muons that are provided by the 
accelerator beam line. The dataset described in this work 
was obtained during the /x^He"'’ measurement campaign 
in 2013 that offers multiple transitions in the low keV 
x-ray region. These consist of the Lq,, and L.y tran¬ 
sitions at 1.52 keV, 2.05 keV and 2.30 keV respectively 
as well as the Kq,, and K.^ transitions at 8.22keV, 
9.74 keV and 10.28 keV, emitted by the muonic helium 
ions during the so-called atomic cascade within a time 
frame of few ns total^^. 

The muons decay after an average lifetime of 2.2 /xs into 
muon neutrino, electron antineutrino and "high energy" 
electrons in the MeV range. These electrons deposit en¬ 
ergy when transversing the APD, creating electron hole 
pairs in all regions of the APD quasi simultaneously. The 
induced signals correspond to virtual x-ray energies of 
up to 50 keV. This would raise background effects for the 
experiment that uses the recorded 8.2 keV ^^He^ Kq 
x-rays as signal for laser spectroscopy. Therefore a sup¬ 
plemental set of 4 plastic scintillators surrounds the gas 
target and APD arrays for additional means of electron 


detection and exclusion of background. Since the overall 
detection efficiency for electrons in the mentioned plastic 
scintillators is only roughly 30 %, additional means for 
electron identihcation were desirable. This was achieved 
by waveform analysis described in Sections III and V. 
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